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0\ . Short -range columnar defects essentially influence the magnetic properties of a mesoscopic su- 

perconducting disc. They help the penetration of vortices into the sample, thereby decrease the 
sample magnetization and reduce its upper critical field. Even the presence of weak defects split a 
' . giant vortex state (usually appearing in a clean disc in the vicinity of the transition to a normal 

P5 ' state) into a number of vortices with smaller topological charges. In a disc with a sufficient number 

, of strong enough defects vortices are always placed onto defects. The presence of defects lead to 

' ' ■ the appearance of additional magnetization jumps related to the redistribution of vortices which are 

0^ ' already present on the defects and not to the penetration of new vortices. 

^ I 

■ PACS: 74.60. Ge; 74.60.Ec; 74.62.Dh 

■ 

O ■ I. INTRODUCTION 

U • 

Advances in microtechnology have allowed the fabrication of Hall probes of micronpgi^c- They were successfullM 
^ ■ applied for time- and space-resolved detection of individual vortices in superconductoril'Bflu. Recently Geim et. al.Q 
^ X ' developed Hall probe techniques by employing submicron ballistic probes of this type for studying individual submi- 
, cron samples. The use of Hall probes in the regime of ballistic electron transport and samples of size smaller than the 
' probe size allowed them to make a link between the detected signal and the sample magnetization. The experiments 

1 showed that the sample undergoes a sequence of phase transitions of the first kind, which manifeststhemselves by 
■ O mesoscopic jumps of the magnetization curveO. These jumps are due to penetrations of additional vortices inside the 
^ superconductor as the applied magnetic field increases. (Due to the small size of the sample, each vortex carries a 

magnetic flux smaller than a single superconducting flux quantum $o-) 



o 



The results obtained in Ref.i stimulated a series of theoretical worksiiiSQ. Deo et. al. numerically solved 
' the 3D non- linear Ginzburg-Landau equations together with the Maxwell equations. They emphasized role of finite 
' sample thickness and showed that S-N transition in mesoscopic disc could be first or second order. They also anal- 
' ysed the conditions of multi-vortex.states or a. giant vortex state formation, constructed a vortex phase diagram and 
\ expalined the experimental resultsD. PalaciosEI considered the same problem within a variational approach, obtained 
T-H ■ the magnetization jumps related to the penetration of new vortices into the sample and showed that below the upper 
critical field for an infinite sample Hc2 the vortices occupy spatially separated positions (a vortex glass structure) while 
, above this field they always form a giant vortex located at the-disc center. Choice of the GinzbuEgrLandau parameter 
' K — i resulted in a good agreement with experimental resultgj. Recently Akkermans and MallikEJ considered a finite 
sample at the dual point with Ginzburg-Landau parameter k — l/\/2 and obtained the magnctizatinH curve which 
also came to the qualitative agreement with the numerical results of RefJ3 and the experimental curveEl. 



. Introduction of strong pinning centers such as columnar defects, which can be produced by heavy-ion irradiationllj, 

(~| ' essentially influuences the magnetic prcmerties of the sample. In bulk superconductors these defects lead to important 

O 1 change of the reversible magnetizationEj. Even small concentration of defects modifies the magnetization curve of a 

J-^ ■ conventional superconductor near Hc2 leading to a sequence of reentering transitions related to the two possible types 

^ , of the local symmetry near each defecto. 
•I— I , 

Columnar defects should also essentially change the magnetic properties of mesoscopic superconductors. Such de- 
fects are known to be insulating inhomogeneities. Generally they can be described as local inclusions with lower 
■ ■ ■ critical temperature. In the case when the number of defects is of the order of the number of vortices one can expect 
that they will essentially suppress the magnetic response of the sample and reduce the upper critical field _ffc3 • If the 
number of defects is larger than the number of vortices and the defects are strong enough it seems plausible that all 
vortices could be pinned by defects. As the applied field changes the vortices can change their position on the defects. 
These rearrangements should lead to increasing of the number of mesoscopic jumps of the magnetization curve as 
compared with that of a clean sample. In the present paper we show that all these phenomena really take place in 
small enough superconducting discs. 
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The content of the paper is as follows. In the second section we formulate the problem. Section NI^ has an auxuliary 
character - here we reproduce some numerical results which should be used later on. In section we describe the 
variational approach for the thermodynamic potential. Properties of the clean disc are discussed in sec tion ^ The 
main results concerning the disc with defects are presented in section VI and summarized in section VII . 



For convenience any length appearing below is measured in units of the temperature dependent coherence length 
^(T). In these units the penetration length coincides with the Ginzburg-Landau parameter n. 



II. THE MODEL. 



Consider a type II superconducting disc with thickness d and radius tq containing columnar defects of size I. The 
sample is subject to an applied magnetic field, which is parallel both to the defects and to the disc axis. In what follows 
we use the dimensionless variables measuring magnetic field and vector potential in units of Hc2 = *J'o/27r^^(T), and 
^q/2tt^(T) respectively ($o is the superconducting flux quantum). Then the density of the thermodynamic potential 
and the order parameter will be measured in units ot^/ (i and v^— ao//3 where ao and (3 are the standard Ginzburg- 
Landau coefficients of the clean disc. In the presence of defects the coefficient a should be modified and depends on 
coordinates 



a(r) — ao(l ^ (5a(r)). 

The last term in parenthesis is simply related to the critical temperature change 5Tc{t) caused by defects: 
where Tc is the critical temperature of a clean sample. 



We assume that the disc is thin and small d <C tq < k. All the dimensions of such a disc are smaller than the 
penetration depth k. Therefore the problem becomes essentially 2D one, and, moreover, it is possiblej-peglect the 
spatial variation of the magnetic induction b inside the disc and replace it by its average value (b)Ll (here and 
below the brackets (..) mean averaging over the sample area). As a result one gets the following expression for the 
thermodynamic potential density 

G = (-m^ + il*!'* + |D_*|2 + <5a(r)|^|2^ + K'{{b) - hf. (2) 

The gauge invariant gradient D_ is given by 

D = i— +a, 
or 

where the symmetric gauge a = {b)r/2 is adopted. 



According to the general approach of the Ginzburg-Landau theory one has to minimize the thermodynamic potential 
density (H) with respect to the order parameter ^E" with an average induction (6) fixed and then to minimize the result 
once more with respect to (6). The first step results in a nonlinear differential equation with a boundary condition 

D_ = 0, (3) 

the solution of which is rather difficult even in the absence of defects. Therefore we use the variational procedure 
choosing the trial function as a linear combination of the eigenfunctions of the operator (D_)^ with the boundary 
condition (||). The corresponding eigenfunctions A^^^ and eigenvalues (T„,m depend on the disc radius rg. Here m is 
an orbital number and n stands for the number of the Landau level which this eigenvalue belongs to when the disc 
radius vq tends to infinity. In strong enough magnetic field one can take into account only n = states and therefore 
the quantum number n will be omitted in what follows. For an infinite sample such an approximation corresponds to 
neglection of higher Landau levels contribution which is justified from the fields h = 0.5tJ. In our case it is adequate 
when the strength of defects Sa{r) is much smaller than the distance between the n = and n = 1 eigenvalues. Then, 
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to describe states with a fixed number A^^, of vortices the maximal orbital number or topological charge which enters 
the trial function should be equal to N^. Finally our trial function can be written as 



* = H exp(-imi9) A„, (4) 

m=0 

where Am is given by 

Am = yWexp (-yW) $ (^^^^'"^+ l^^^O ■ 

In Eq.(^) the expapsion coefficients Cm serve as variation parameters and <i>(a, c; x) in Eq.(^ is the confluent hyper- 
geometric functionllB. 

To proceed the problem one should substitute the trial function into the expression for the thermodynamic 
potential density and first minimize it with a respect to the expansion coefhcients Cm at an average induction (6) 
fixed. As a result one obtains a system of a finitfi number of nonlinear equations for the coefficients Cm- This system 
is a finite version of the Ovchinnikov equationsU. However in the presence of disordered set of defects the solution of 
thcsp-|eq|Wations is very complicated. The point is that now no selection rule (successfully used in the homogeneous 
caseE-lllJa) can be applied. Thus the problem needs another approach. 

In what follows we consider a disc which contains Nd short-range defects I <C 1 placed at the points ri, r2, tat^. 
The number of defects is assumed to be larger than the maximal possible number of vortices . As we could see 
(see section ^ below) a small enough clean disc can accumulate vortices only in its center. The defects attract the 
vortices and due to their short range can pin the latters exactly on the positions of the defects. Therefore we consider 
only some special configurations of vortices such that they occupy only the positions of defects and the disc center. 
This choice of trial function implies the following procedure. Let us fix a defect configuration {rj}, j = 0, 1, ...,iVd, 
ro = 0, a set of corresponding topological charges {p(j)}, an external magnetic field h and an average induction (b). 
Each topological charge rij is non negative integer and the set {p{j)} satisfies the condition 

^Pij) = N.- (6) 

Thus our procedure accounts for the existence of multiple vortices located on the disc center or on any defect position 
as well. The trial function (^ has zeros only at points {rj} with miltiplicities p{j). The latter condition completely 
defines all coefficients {Cm} i'm- — 0, 1, -/V^, — 1) up to a common multiplier Cn^, which we term as the order 
parameter amplitude. Further, we need to minimize the thermodynamic potential with respect to this amplitude and 
the average induction. The result has to be compared with those obtained for different total numbers of vortices and 
different sets of "occupation numbers" {p{j)}- Comparing the obtained value of the thermodynamic potential with 
that corresponding to a normal state one finally finds the preferable state of the disc for a fixed value of external mag- 
netic field. Repeating this procedure for various values of the magnetic field one could describe magnetic properties 
of the sample in a wide range of the fields up to the upper critical field Hcs ■ The next four sections are devoted to 
the realization of thcd procedure described above and to the presentation of its results. 



III. SPECTRUM OF THE OPERATOR (D^)^ 

To construct the trial function (^ one should first obtain the eigenvalues and eigenfunctions of the operator (D_)^. 
This is a textbook problem and it was solved many times but we need the solution for various disc radii and various 
average induction values. The eigenvalue equation reads: 



1 d ( 1 / d 



VT — 
dr } ' r'^ Kd-Q 2 



2 

\h)rA * = (7) 



dr 



(8) 
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Solution of this differential equation can be written as: 



where C = (6)r^/2. The boundary condition (|^) implies the following eigenvalue equation for the quantities cTn^m' 



Co V 2(') " ' 7 (4)(l"l + l) V 2(1,) 

where C,q — (&)rQ/2 and the index n stands for the number of a Landau level, to which the quantity (Jn.m tends as the 
disc radius tends to infinity: 

lim Gn.m — {b){2n + \m\ — rn + 1). 

ro — >oo 

We solved equation numerically tabulating some needed eigenvalues cr„,m and the corresponding eigenfunctions 
Am for various quantum numbers n = 0, 1, m = 1,2, 3, 4, 5 and disc radius r = 2.6. The eigenvalues as functions of an 
aver|age induction are shown in the fig. |l]. These results are completely consistent with e.g. those obtained earlier in 
Ref.E£l. One can observe that the distance between the zeroth and the first Landau levels is of the order of unity. So 
we can indeed neglect in expansion (^) the contributions of higher "Landau levels" as long as defects are not extremely 
strong, 6a{r) < 1. 

The results shown in fig. |l| help us estimate how many vortices can enter the sample. Indeed, for a = 1 the eigen- 
value equation (|^) coincides with the linearized Ginzburg-Landau equation. Therefore the maximal average induction 
{b)m corresponding to (7™ = 1 can be treated as the upper critical field for a given orbital number m. The highest 
of these fields is the genuine upper critical field /ic3 and the corresponding value of m gives the topological charge of 
the giant vortex usually appearing in the vicinity of the clean disc phase transition point (see Refs .11313 and section 
below). In the case r = 2.6 considered here the highest possible field at which superconductivity still exists is 
hc3 ~ 1.98. This corresponds to the intersection point of the curve (74 and the dashed line a — 1. Thus a clean 
superconducting disc of this radius at the phase transition point can accumulate only four vortices since the curve for 
n = 0, m = 5 never reaches the line a — 1. 




FIG. 1. Eigenvalues an,m for the disc of radius ro — 2.6 as a function of the applied field h. 



IV. THE THERMODYNAMIC POTENTIAL 



Substituting the test function for the order parameter into the expression for the thermodynamic potential 
density (0) one obtains: 



G = - |C™P (1 - a™) lm + ^ C:,,C:CkC^+n-kJm,nM + {Sal'fl^) + - hf 



(10) 



A;,m,n— 



where the brackets (..) mean averaging over the sample area, Im = (^m)' Jm,n,k = (AmA„ Am+„_fc) and am = 
co.m- For the state characterized by a topological charge Ny the coefficient Cn„ necessarily differs from zero. We choose 
it as an amplitude of the order parameter and introduce new expansion coefhcients Dm and new order parameter ^ 

Cm — Cpf^Dmj 

Dn^ = 1, (11) 
= Cn^iP- 

Rewriting the thermodynamic potential (|l^) in terms of these new variables and varying it with respect to the 
amplitude Cn^ we obtain the following expression for its extremal value: 

^(1-a™)/™- (fc|7/-|2) 
\CnJ'^ =f2 . (12) 



(6) ^ D*^D*^DkDm+n-kJn 



fc,m,n— 



The expansion coefficients of the order parameter (^, (|11|) are completely defined by the position of vortices on the 
defects. Let us choose some configuration of vortices {r^}. In this set there are points occupied by a single vortex 
{p{j) — 1) and points corresponding to multiple vortices with topological charge p{j) > 1. Then the set of coefficients 
{Dm} — {f~^Cm} can be calculated from the following system of linear equations: 



, (p(i)) , 



J2 i?™exp(-jm7?,)A(PW)(r,) = exp(-ziV,z9,)Al^y" (r^), (13) 

where the notation /'-"'(x) is used for the nth derivative. 

The "inhomogeneous term" in Eq.(|^) which is proportional to 5a{r) appears due to columnar defects. We have 
already mentioned that defects are supposed to be short-range ones. In this case this term can be represented as a 
sum over defects. For the Gaussian form of defects 

.aW^5ii:c.xp(-lL_ii)!) 

j=l \ / 

the "inhomogeneous" term in (|lO|) in the leading approximation with respect to our small parameter I can be rewritten 

as 

2 

Substituting equations ( pi] ) , ( p2| ) and (|l^) into Eq. ( |lO| ) we obtain the final expression for the thermodynamic potential 
of the disc with defects: 

\ 2 

Nd 

Y^\Dm\'{l-'Jm)In.-^Y.\^{r,)f' 

G = + n\{h) - hy. (16) 

2(6) ^ DmD*^DkDm+n~kJm,n,k 
k,m.n—0 

We solve the system ([l3| ) for each combination of vortices on the defects in order to find the set of expansion 
coefficients {Dm} as a function of the average induction (b). The set of coefficients is then pluged into expression ([l^ ) 
for the thermodynamic potential G at a fixed applied field h. Now we can find the average magnetic induction (b) 
at which the thermodynamic potential (^^ has a minimal value at fixed applied field and configuration of vortices. 
After that we must repeat this procedure for different configurations and different values of the applied field. As a 
result, we obtain a number of data sets for the thermodynamic potential as a function of the applied field for different 
configuration of vortices. Then for each value of an applied field we should choose the preferable vortex configuration 
which minimizes the thermodynamic potential. This enables us to obtain the disc magnetization as a function of the 
applied magnetic field. 
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V. CLEAN DISC 



We start from the case of a clean disc with radius vq — 2.6 and k = 3. Ahhough this value of k limits the condition 
K'3> To, the chosen region of applied fields enables us to neglect the spatial variation of the magnetic inductionB.The 
maximal number of vortices in such a disc equals four (see section HI). Due to the sample geometry and small 
maximal number of vortices they can form only a number of symmetric configurations when some vortices occupy 
the disc center and the others are placed away from the center in such a way that they form a regular polygon. All 
these configurations are presented in fig. 2. In cases (b), (h); (d) and (g) the topological charge of the multiple vortex 
at the origin is equal to 2; 3, 4 respectively. In cases (c),(e),(f),(h),(i),(j) the shifted vortices are place at a distance p 
from the origin. 




:.)0(3 




FIG. 2. Possible configurations of vortices inside clean disc of radius — 2.6. 

For a given vortex configuration the expansion coefficients {Dm} can be calculated from the system of linear 
equations (^. For each possible vortex configuration we substitute these coefficients into the expression for the 
thermodynamic potential of the clean disc 



iA„i'(i-fT„o/„, 

,m=0 / 



+ K\{b)-h)' 



2{b) D*mDlDkDm+n-kJrr 



k,m.n—0 



and minimize it with respect to the average induction (6). We repeat this procedure for all configurations and for 
various distances of vortices from the disc center inside each configuration. Thus the problem has three variational 
parameters: the type of vortex configuration (fig. 2), the distance p of vortices from the disc center and the average 
induction (6). We changed the distance p by step of 6p = O.lro- Numerical calculation showed that because of the 
disc small size only configurations in which p = (fig. 2 (a),(b),(d),(g)) gain the energy. So within the calculation 
accuracy 6p — 0.26 we have only a multiple vortex at the disc center with a possible topological charge p(0) = 1, 2, 3, 4. 
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FIG. 3. Magnetization curve of a clean superconducting disc of radius ro = 2.6. 

The dimensionless magnetization m ^ h~ {b) of the clean disc is presented in fig. ||. Penetration of an additional 
vortex inside the sample is manifested by magnetization jump. Each branch of the|4:urve corresponds to the one-, 
two-, three- and four- vortex states. This result is similar to that obtained by PalaciosO and Deo et. alu for discs with 
larger radii and it will be used in th next section devoted to the magnetic properties of the disc with defects. 



VI. DISC WITH DEFECTS 

In the case of disc with defects, one should take into account the defects configuration and minimize the thermody- 
namic potential . We present below the results for a single configuration of the defects obtained with the help of 
a random number generator. We hope that it is rather typical (see fig. ||). In any case the results obtained below for 
this configuration enable us to demonstrate all the new features characterizing the magnetic properties of a sample 
with defects and to confirm all the expectations formulated above in the Introduction. 




-2.6 -2.0 -1.3 -0.6 0.0 0.6 1.3 1.9 2.6 

FIG. 4. Defects positions in the disc 
The coordinates of defects are collected in Table |. (Note that all distances are measured in the temperature de- 
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pendent coherehce length units.) 



We analyze the thermodynamic properties of the disc for various values of defect strength ai. This constant can 
be easily varied experimentally by changing the sample temperature (see Eq.(0)). To present the results more clearly 
we collect all configurations of vortices which will be realized for values considered for the defect strength in Table 0. 

The left column of the table contains the values of the coupling constants. The upper line enumerates the vortex 
configurations ordered with accordance to their appearance with the growth of a magnetic field. The same numbers 
enumerate different regions of the magnetization curves on figs. ||, ^ Note that the last configuration in each line 
appears just before the phase transition to the normal state at the upper critical field hcs- Then, each configuration 
is described by an ordered sequence of six numbers. The j-th number is equal to the topological charge located at 
the point Tj-i. In other words the first number is the topological charge at the disc center, the second number is the 
topological charge at the first defect and so on. For example configuration {211000} corresponds to double vortex at 
the disc center and two single vortices placed at the first and the second defects. 



TABLE I. Coordinates of Defects. 





X 


y 


r 


1 


-0.253 


-1.755 


1.773 


2 


0.830 


0.856 


1.192 


3 


-1.205 


-1.248 


1.734 


4 


-0.755 


0.948 


1.212 


5 


1.083 


1.405 


1.774 



TABLE IL Configurations of vortices. 





1 


2 


3 


4 


5 


6 


0.04 


100000 


200000 


300000 


40000 






0.08 


100000 


200000 


300000 


310000 






0.12 


100000 


200000 


300000 


211000 






0.16 


100000 


200000 


101100 


300000 


101200 


211000 


0.3 


000101 


000110 


000111 


001110 


000130 


001210 
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We start from small values of the defect strength. The corresponding magnetization curves are shown in fig. 




FIG. 5. Magnetization curve of the superconducting disc of radius ro — 2.6 and k = 3 in the presence of defects with effective 
coupling constants ai — 0.04 (a), qi — 0.08 (b), ai — 0.12 (c) and qi = 0.16 (d). 

The first part (a) of this figure describes the magnetization curve for a sample with ai — 0.04. Because of the 
small value of the coupling constant, this part is qualitatively equivalent to that for a clean disc. Each branch of 
the magnetization curve corresponds to a one-, two- , three- and four-vortex states. These branches are divided 
by jumps of the magnetization which are caused by penetration of an additional vortex inside the sample. However, 
even in this case some new features caused by defects are manifested. We particularly refer to the suppression of 
magnetization, penetration of new vortices at lower fields and decreasing of the upper critical field in comparison 
with the results for the clean sample (see fig. ^). Magnetization of the samples with ai = 0.08 (fig. ^.b) and with 
tti = 0.12 (fig. I^.c) have the same number of mesoscopic jumps as in the previous case. This means that all the jumps 
are still due to vortex penetrations. However a new interesting feature appears near the phase transition point. The 
four-multiple vortex at the disc center is split. In the case ai = 0.08 (fig. ^.b.4) three-multiple vortex remains at the 
center and one more vortex occupies the first defect (configuration {310000}). The corresponding distrubution of the 
absolute value square of order parameter is presented in fig. |l 




FIG. 6. Square modulus of the order parameter for ai = 0.08 at an applied field h — 1.753. The vortex configuration is 
{310000}. 

More complicated splitting is observed in the case ai = 0.12 (fig. ||.c.4) Two vortices remain at the disc center, one 
occupies the first defect and another one occupies the second defect (configuration {211000}). The square modulus 
of the order parameter is plotted in fig. ^. 



9 



FIG. 7. Square modulus of the order parameter for ai — 0.12 at an applied field h = 1.7. The vortex configuration is 
{211000}. 

In the two latter cases the defect strength was relatively small. Therefore the defects could partially destroy the 
giant vortex state with maximal multiplicity which precedes the transition to the normal state. Further increasing of 
the coupling constant leads to appearance of additional mesoscopic jumps related to the rearrangement of the vortices 
on the defects as the applied magnetic field changes. Consider the case ai = 0.16 (fig. ^.d). At small values of 
the applied field one gets one- and two-vortex states at the disc center. However, when the third vortex is allowed 
to penetrate (fig. |[d.3) the multiple vortex is destroyed and the vortices occupy the disc center, the second defect 
and the third defect (configuration jlOllOO}). Plot of the square modulus of the order parameter for this vortex 
configuration can be found in fig. H (to present the plot more clearly the orientation of the axes is changed with 
respect to the two previous plots). 

With further increasing of the applied field the system turns again into the three-multiple vortex state at the disc 
center (fig. ^.d.4). So in the same sample two different vortex configurations with the same total topological charge 
are possible. When the fourth vortex penetrates the disc the three-multiple vortex state splits again (fig. |.d.5) into 
double vortex at the third defect, one vortex at the disc center and another one at the second defect (configuration 
{101200}). The appearance of the second vortex on the third defect is a result of a very restricted space of the trial 
functions. Indeed, according to Eq. (^6|) any defect which is already occupied by a vortex is put out of the game and 
one can not gain energy adding one more vortex to the same defect. This means that in a wider variational space the 
configuration {101200} would be replaced by another one which should be more preferable. At the same time it will 
necessary lead to the corresponding magnetization jump. 

With increasing of the applied field we have a new jump of the magnetization curve, which is caused by rearrange- 
ment of the vortices into the configuration {211000} identical to that of the four vortex state in the case ai — 0.12. 




FIG. 8. Square modulus of the order parameter for a^ff — 0.16 at an applied field h = 1.4 The vortex configuration is 
{101100}. 

Thus one can see that the stronger defects are the greater is the tendency of vortices to occupy defects. The 
destruction of the giant vortex at the disc center begins near the upper critical field. Increasing the defect strength 
destroys the centered multiple vortices with lower multiplicity. The preferable arrangement of the vortices corresponds 
to the maximal reduction of the square order parameter modulus. 

At strong coupling constant one expects to get states where all vortices are placed onto defects for all values of the 
applied field. Consider the results of studying the case ai = 0.3. The magnetization curve of such disc is shown in 
fig. |. 
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FIG. 9. Magnetization curve of the superconducting disc of radius tq = 2.6 and k = 3 in the presence of defects with an 
effective coupling constant ai =0.3. 

Penetration of vortices inside the disc with such strong defects occurs at values of the applied field smaller than 
that of the previously considered discs with relatively weak defects. Because of that, already at a field h = 0.6 the 
disc accumulates two vortices (fig. ^.a,b.l). Their configuration is {000101} (see fig. p^ . 



FIG. 10. Square modulus of the order parameter for ai — 0.3 at an applied field h = 0.71. The vortex configuration is 
{000101}. 

As the applied field increases this configuration is changed by another one {000110} with the same total topological 
charge. Three vortices appearing at higher fields always occupy three different defects. The corresponding configura- 
tions are {000111} and {001110}. Two configurations with total topological charge four are realized. Both contain a 
multiple vortex on one of the defects. The first configuration appearing in relatively low field is {000130}. Here one 
has three-multiple vortex on the fourth defect. The second configurarion {001210} preceding the transition to the 
normal state at hc3 contains a double vortex at the third defect. Plots of the square modulus of the order patameter 
for these cases are shown in figs. |ll| and 12. Thus in the case of a strong defect ai = 0.3 considered here the number 



of magnetization jumps within the same field region is twice the number of possible values of the total topological 
charge. We do believe that in a disc of the same radius containing more defects this number will increase. 



FIG. 11. Square modulus of the order parameter for ai — 0.3 at an applied field h = 1.31. The vortex configuration is 
{000130}. 
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FIG. 12. Square modulus of the order parameter for ai — 
{001210}. 



0.3 at an applied field h = 1.55. The vortex configuration is 



We already mentioned that the presence of attractive defects reduces the upper critical field hcs at which the 
thermodynamic potential of the superconductor ( |l6|) becomes equal to zero (the thermodynamic potential of normal 
metal) . Figures ^ and |9| show that the larger the defect strength ai is the lower is the transition field. The dependence 
on the upper critical field of the defect strength ai is shown in fig. |l^. 



o o 



o o 



FIG. 13. The upper critical field as a function of the defect strength. 



VII. SUMMARY 



Summarizing, we studied magnetic properties of mesoscopic superconducting discs with disordered attractive colum- 
nar defects. The number of defects is assumed to be larger than the maximal possible number of vortices accumulated 
by the disc. We obtained the magnetization curves for various strengths of defects in a wide region of the applied 
magnetic field. The results show that the defects help the penetration of vortices into the sample. They also reduce 
both the value of the magnetization and the upper critical field. Even the presence of weak defects can split the giant 
vortex state at the disc center (usually existing in a clean disc of small radius) into vortices with smaller topological 
charges. This splitting occurs in the vicinity of the upper critical field. Strong chough defects always pin all vor- 
tices, splitting multiple vortex states at the disc center in all field region. This leads to the appearance of additional 
mesoscopic jumps in the magnetization curve related not to the penetration of new vortices into the sample but to 
redistribution of vortices within the set of defects. The number of these jumps enlarges increases with the number of 
defects. 
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